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In this paper we derive upper and lower bounds of the spectrum of a 
sinusoidal carrier phase modulated by gaussian noise having a rectangular 
power spectrum. It has been jound in practice that such a random process 
adequately simulates for some purposes, a frequency division multiplex 
signal, a composite speech signal, and so on. We show that these upper and 
lower bounds of the spectrum are very close to each other if the root mean 
square phase deviation of the carrier is even moderately high. Also, a simple 
method called the saddle-point method can be used at all frequencies f to 
estimate the spectrum with less than ten percent error. We also show that 
the results obtained from the quasistatic approximation, often used in 
such cases, are too small for large f, and that this low-frequency approxi- 
mation cannot be used in cases where the behavior on the tails is 
of importance. 

I. INTRODUCTION 

It has been found in practice that a bandlimited random gaussian 
noise having a rectangular power spectrum adequately simulates for 
some purposes a wideband composite speech signal, a frequency divi- 
sion multiplex baseband signal consisting of a group of single side- 
band carrier telephone channels, and so on. 1 In the design of commu- 
nication systems, the spectral characteristics of a sinusoidal carrier 
phase modulated by such a baseband signal are of great interest; 
various methods have been used in recent years to estimate this 
spectrum for large and small values of mean square phase deviation 
of the wave, both close to and far from the carrier frequency (that is, 
in the principal part of the spectrum and far down on the tails of the 
spectrum respectively) . 1_8 

It has been shown that the spectrum may be expanded as an 
infinite series of weighted convolution terms. 2 ' 5 - 7 ' 8 This series may 
be used to estimate the principal part of the spectrum (close to the 
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carrier) for small or moderate index (that is, small or moderate values 
of rms phase deviation). However, for large index, or far down on 
the tails for small index, too many terms would have to be included 
if this series is to be used directly. 

The simplest analysis — often called the quasistatic approxima- 
tion — yields a gaussian spectrum for large-index angle-modulated 
waves 2, *- 8 in most cases.* This approximation fails far out on the 
tails of the spectrum; a careful investigation has been given in only 
a few cases. 7 We obtain below upper and lower bounds for the spec- 
trum of an angle-modulated wave with white, band-limited phase 
modulation; far out on the tails the spectrum far exceeds that pre- 
dicted by the quasistatic approximation. 

This problem is of interest in considering interference between 
two (or more) phase modulation (PM) systems in neighboring loca- 
tions. Consider the following situation. In the frequency bands above 
10 GHz, where the signal attenuation due to rain storms could be very 
severe, close spacings of the repeaters are almost mandatory for 
reliable communication from point to point. In such cases the prob- 
lem of interference between neighboring systems may be much more 
important than the problem of noise; the system may thus be in- 
terference limited. In order to combat this interference it is very 
likely that broadband modulation techniques like PM [or frequency 
modulation (FM) ] or pulse code modulation (PCM) will have to be 
used. In order to investigate the effect of this interference between two 
co-channel PM (or FM) waves it is necessary to evaluate the spec- 
trum of a PM wave, so that the parameters (such as rms phase 
deviations) of the two PM systems can be properly chosen to keep 
the interchannel interference below a certain desired level. 

We first obtain an expression for the covariance function of the 
PM wave. From this covariance function we then derive an expres- 
sion for the spectrum of the PM wave and show that it can be 
expressed as an infinite series. This series has been evaluated for 
certain values of rms phase deviation iV. 6 

We then show that the autocorrelation function of the PM wave 
is analytic at all points in the finite part of the complex plane deter- 
mined by the argument of the autocorrelation function. In determin- 
ing the Fourier transform of the autocorrelation function we change 
the path of integration! so that the contour is very close to the path 

* For exceptional cases see Ref. 7, Ch. 4, pp. 131-135. 

t The method used in this paper to evaluate the spectrum is a close relative 
of the method of steepest descent (or the saddle-point method) used in evaluat- 
ing certain kinds of integrals. 1, 7 ' 10, u - M 
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of steepest descent of the integrand. We then divide this contour into 
four (or five) disjoint sections and show that the major contribution 
to the integral comes from one of these sections. 

We next derive upper and lower bounds to the spectrum S v (j) of the 
PM wave and show that these bounds are very close for all / and for all 
values of rms phase deviation N ^ 5. For N ^ (10)* we show that the 
spectrum can be evaluated by this saddle-point method in a very 
simple manner with a very small fractional error (less than 10 percent), 
and we give this method. 

We finally compare the quasistatic approximation to the saddle- 
point approximation. For large values of frequency /, we show that 
the quasistatic approximation gives too small a value for the spec- 
trum, and that it cannot be used in cases where the spectral behavior 
on the tails is of importance. However, as we show, the saddle-point 
method can be used in all cases in which N is moderately high. 

In conclusion, this paper gives a simple method of evaluating the 
spectrum of a sinusoidal carrier phase modulated by gaussian noise 
having a rectangular power spectrum and having a moderately high 
modulation index. 

II. SPECTRAL ANALYSIS OF PM WAVES WITH RANDOM PHASE MODULATION 

A sinusoidal wave of constant-amplitude phase modulated by a 
signal n(t) may be written as 

W{t) = A cos [wo< + n(t) + 0], (1) 

where A is the amplitude of the wave, / = o) /2ir is the carrier fre- 
quency of the wave, and 6 is a random variable with probability density* 



TTo(6) = . 



i ■ ° s • < ^ ( 2) 



0, otherwise. 

Assume that n(t) is a stationary bandlimited white gaussian ran- 
dom process with mean zero and variance N 2 .~\ Its spectral density 



*If n(t) is a stationary random process the introduction of random variable 
in equation (1) makes W(t) a random process which is at least wide-sense 
stationary so that its spectrum can he calculated from the Wiener-Khintchine 
theorem. 2,7 If we do not have 6 in equation (1), W(t) is no longer (even wide- 
sense) stationary, and the spectrum of Wit) is usually calculated from the time 
autocorrelation function of W(t)? The results obtained in the two cases are 
identical. 

t The parameter N represents the rms phase deviation (or modulation index) 
of the PM wave given in equation (1). 



Sn(f) - 



(3) 
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S„(/) can be represented (see Fig. 1) as 

[N 2 /2W, I / I < W, 

0, otherwise. 

Such a random process n(t), is often used to simulate a multiplex 
signal, a composite speech signal, and so on. 1, B 

We can show from equation (3) that the covariance function R n (r) 
of n(t) is given by 



Rn(r) = N' 



sin 2ttWi 



(4) 



2tWt ' 

this function R„(t) is shown in Fig. 2. Since n(t) is a stationary gaus- 
sian random process it can be shown that W(t) is at least wide-sense 
stationary and that its covariance function R w {t) can be represented 
as 2 - 7 



A 2 
R w {t) = — exp [—R N (0)] exp [R n (t)] cosw t. 



(5) 



From the Wiener-Khintchine theorem, and from equation (5), the 
spectrum S w (/) oiW(t) can be written as 



or 



where 



S W (D = [ Rw(t) exp [-;2tt/t] dr, 

J —Ml 



Sw(j) = ~ [Sy(f ~ /o) + S v (f + /o)], 



S v (f) = f exp [-R N (0)] exp [R„(r)] exp \-j2tJt) dr. 

J -to 



(6) 



(7) 



(8) 



SPECTRAL DENSITY 
IN RAD 2 / HZ 



2W 



-W W FREQUENCY IN Hz 

Fig. 1 — Spectral density of phase modulation. 
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Fig. 2 — Covariance function of nit). Since R n (r) is an even function of r we 
only show R n (r) for t 5^ 0. 



From equations (4) and (8) we have 



^-sarL^r 



■A^ s 



1 ^J| exp [-jXp] dp, (9) 



.sin p 



where 



X = 



TF 



(10) 



III. SERIES EXPANSION OF SPECTRUM FOR GAUSSIAN MODULATION 

The integral in Equation (9) can be evaluated by expanding 

-H>-T]} 

into a Taylor series ; integrating term by term we can write" 

^(4-Tll - - [ -** ! KM'- (,i) 

* We note that £™_o x"/n\ converges uniformly to exp [x] for all finite values of x. 
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From equations (9) and (11) we have* 
S v (i) = exp [-JV 2 ] 

■fa + drSStW -<-**}* (12) 

where 8(f) is the Dirac delta (unit impulse) function. 
We now note that 



/; 



^ exp [-jXp] dp = *\(X) = f' I X ! < *' (13) 

^ lO, otherwise, 

or 

^i(X) = ^[w-i(a + 1) - «-i(* - 1)]. (14) 

where it-i (a;) is the unit step function defined by 

w_,Ct) = \ (15) 



10, x < 0, 
and that 1 - " 

f (*^J exp [-jkp] dp = F t (\) 



P 

^f(-i)' ( ' x ]+/rff- . os|x|<f«*a, (16) 

0, otherwise, 

where 

ilf = JOT ^ + J X ! ] , (17) 

and 7JVT[x] represents the largest integer not greater than x. 

It can be shown that F«(X), / ^ 2 is a continuous function of X and 
that Fi(X) is discontinuous at X = 1. For large t, we can show from the 
central-limit theorem f that 2 

^-(fy^of]- (>8) 

* The term containing 8(f) in equation (12) represents the dc component of S v (f). 

fain p] t 
t See pp. 362-366 of Ref. 2. It can be shown that I — ^ — J can be interpreted 

as the characteristic function of the sum Q of / independent random variables with 
identical uniform probability distributions. 6 The function F t (\)/2ir therefore 
represents the probability density function of fi. Alternatively F t (\) is the (I — 1)- 
fold convolution of the flat spectrum with itself. 
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From equations (12), (13), and (16) we can write 

Sr(f) = exp [-JV 2 ][a(/) +^Z^ F,(//TF)] • (19) 

For N 2 = 6, we have calculated the spectrum from equation (19) 
and the results are shown in Fig. 3. Notice that the spectral density is 
discontinuous at j/W = 1. 

For JV 2 <5C 1 (low-index modulation), we have from (19) 



Sv(f) « S vo (ft = exp [-JV 2 ][a(/) + ^jy F M w )\ • 



(20) 



and the error in this approximation may be investigated from equations 
(9) and (20).* The approximation given in equation (20) represents the 
low-index approximation for the spectrum; this result has been obtained 
by many authors. 3-6 

The series given in equation (19) may be used to estimate the princi- 
pal part of the spectrum (close to the carrier) for small or moderate 
index, since only a small number of terms need to be included in the 
partial sum. However, for large N 2 , or far down on the tails of the 
spectrum for small N 2 , too many terms would have to be included to 
estimate the spectral density. In fact for N 2 ^> 1, or for j/W >>> 1, 
the degree of complexity required in estimating S v (j) from equation 
(19) becomes inordinately high. 

When AT 2 ^> 1, and for low frequencies, several authors have given 1 " 7 
the quasistatic approximation f 

S r (» ft, exp <-*•) m + j^ (£f exp [-£ (±)'] ( 21) 

for the spectrum. The question arises whether equation (21) can 
be used for large /. Since Rw{t) is infinitely differentiable there is 
no simple way (known to the authors) of investigating, for large /, 
the error is this approximation. 7 

In the problem of interference between two neighboring channels 
it is necessary to evaluate S v (f) for large / so that the effect of this 
interference can be determined. As we shall show later on in this 
paper equation (21) gives too small values to <Sv(/) for large /; it is 
therefore essential to have a simple and elegant method (different 
from the series method) to evaluate Sy(f) for large / and for large N 2 . 



* At times the low-index approximation for the spectrum is written as exp 
[-JV 2 ]5(/) + (N»/2irW)Fi(f/W). For W << 1, exp [-JV*] « 1. 
t Note that mean square frequency deviation is N 2 W /Z. 
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Fig. 3 — Spectral density of a PM wave for N' = 6. The discrete part of Sv(j) 
for / = is not shown in this figure. Note the discontinuity in the spectrum at 
l/W = 1. 

Readers who might be interested in the final results without wishing 
to work through the detailed analysis, might skip Section IV of this 
paper. 



IV. SPECTRUM EVALUATION BY CONTOUR INTEGRATION 

Let us now consider the integral given in equation (9). Since 8{f) 
and Fi (f/W) are discontinuous functions of /, let us define an integral* 



8(f) = £ exp [-^]{exp [*" 9 -^] - (l + *■ S -^)} 



- (1 •■:- N 2 *-^ 2 )} vn-y/ni/, 



* In Ref. 14 this integral has been studied by Lewin for X = and X = 1. It 
also occurs in several limiting cases in Ref. 1. It is sometimes referred to by 
the name Lewin's integral. 1 
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or 

S(f) = £ exp [-tf']{exp [tf» S -^] - (l + iV ^)} 

■exp(-j\p)dp. (22) 

Since all F,(/), £ ^ 2 are continuous, it can be shown from equations 
(19), and (22) that £(/) is a continuous function of /. 

From equations (9), (19), and (22) we can then write*" 

S r (f) - exp (-#•){«(/) + |£ [?,_,(/ + TF) - «_,(/ - W)] 

+ ^ReS(/). (23) 

Notice from equation (22) that the integration is carried out along 
the real axis, and that for large \\\ (or ] / \/W) , the final factor of the 
integrand exp (—j\p) is a very rapidly oscillating function of p. From 
Refs. 7, 10-13 notice that in such circumstances the method of steep- 
est descent (or saddle-point method), or one of its close relatives, is 
often useful to get an approximate expression for the integral; we 
shall now apply such a method to evaluate S v (f) . 

In applying this method to the evaluation of an integral with a 
real variable of integration, we must first be able to regard the inte- 
gral as a contour integral along the real axis of the complex plane, 
with an analytic integrand. We note that the integrand in equation 
(22) is an analytic function of p, and that it has no singularities in 
the finite part of the complex plane (defined by p) . From Cauchy's 
theorem it therefore follows that the contour of integration can be 
arbitrarily deformed as long as one end is at p = — oo + ;0 and the 
other at p = oo + jO. 11 

In making use of the method of steepest descent the contour must 
be deformed so that the phase of the integrand remains constant (or 
almost so), while the magnitude of the integrand is small except in 
one or more localized regions, where it varies rapidly. This is usually 
accomplished by deforming the contour so that it goes through one 
or more saddle points. In other cases there may be some additional 
constraints on the contour; 7 then only a portion of the path of steep- 

*Re z and Im z denote respectively the real and imaginary parts of complex- 
number z. 
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est descent through a saddle point may be used in rinding the integral, 
and the deformed contour may not actually reach the saddle point. 
In this case the original integral is usually reduced to a virtually real 
integral whose integrand behaves sufficiently simply on the modified 
path of integration so that an approximate evaluation of the integral 
with rigorous (upper and lower) bounds on the error may be obtained. 

Departures from the strict method of steepest descent occur in this 
paper in that approximate paths of steepest descent are chosen. Al- 
though not quite optimum, they are analytically tractable and serve 
to give useful bounds on the integral under consideration. 

Consider equation (22). Since the integrand in equation (22) is 
an analytic function of p, let us assume that p = x + jy is a complex 
variable, and let us deform the contour so as to obtain the path of 
steepest descent. Since the integrand behaves properly on the contour 
for large | p \ , it is clear that the contour of integration can be de- 
formed in quite a general way in the complex p-plane without modify- 
ing the value of the integral. 

From equation (22) it can be shown that the major portion of the 

integrand 

R(p) m exp [-N 2 ] exp \N 2 ^J exp [j\p] (24) 

has a saddle point on the imaginary axis, with the path of steepest 
descent parallel to the real axis at this point. The location p a = jy 8 of 
this saddle point is given by 

coshy, sinh y. X /_ _ (t) . 

~~V. y 2 . " N 2 ~ WW ' {Zb) 

for a given f/N 2 W, equation (25) can be solved numerically to give 
the required y a . We plot y, as a function of f/N 2 W in Fig. 4, and In 
R(jy„) in Fig. 5. 

Let us now deform the contour so that it passes through the point 
Pa = JVa and is parallel to the real axis at this point. From equation 
(22) we then have 



^.^[-^(coeh'l- 8 -^)' 



Re / (26) 



where 
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/ = exp — 



r^sinhy.l 



£M^ sJ ^H 



1 + N'' 



x + ju 
•exp [j\x] dx. 
Rewriting equation (27) 

I = f G(x, y.) dx, 

J -as 

where 
G(x, y.) = exp [-Q R (x, y.)] exp [jQ^x, y.)] 
sinh y 



sin (x + jy.) 
x + jy. J 



— exp —N 



r^sinhy.lf 



y. . 
where Q B (x, y B ) and Qi(x, y„) are real and 

«.(*. v.) - *#£* - Re 



l ^ sin( * + >^ exp(^), 



[" sin (x + J2/.) l 
L z + jy, . 



(27) 



(28) 



(29) 



(30) 
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Fig. 4 — Location of saddle-point j/«. 
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Fig. 5 — Value of -[In R(jy.)V(.2N*). 



or 



1 - 



y. sin x 



and 



tanh y, x , 1 — cos x 

n i \ »r2 sinh y, _ 2 y] x 2 

Q R (x,y.)=N -jf-x l + ( x /y.f 



Q*. V.) = N> Im [*ft^] + ^ 



or 



Q/(z, yO 



sinh ?/, cos x — y, cosh t/. 



= AT 1 



s in a: 



(31) 



(32) 



x- + y, 



+ 



[£2^ i _5^J| JCt (33) 



The functions Q fl (x, y.) and Qj(x, y.) have been plotted in Fig. 6 for 
a set of values of y. . 

Since we are primarily interested in the high-index case let us assume 
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that N 2 » 1. From equations (27)-(33) we now observe that the 
principal contribution to the integral / comes from small x. For small x, 

G(x,y.)ttexp[-Q B (x,y.)]. (34) 

It can be shown (see Figs. 6, 7, and Appendix A) that Q R (x, y.) is a 
monotonically increasing function of x for ^ x fS t and that it oscil- 
lates for values of x > t. For large x, 
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Fig. C — Functions Q R (x, y.) and Qj(x, y,). 
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Fig. 7 — Function Q R (x, y>). From this figure, it can be seen that Q R {x, y.) is a 
monotonically increasing function of x for ^ x ^ it. 



r^sJnhj/TI 



(?(&, y.) « exp 

• s exp iV 2 cosh ?y s 

- 1 + AH cosh y. 



sin g 
a; J 

sin .r 



exp 



jTV 2 sinh y, 



cos a: 



+ j sinh j/. 



COS X 



(35) 



and we note that the first and second terms both have small and 
almost equal magnitude, and almost opposite phase angle, so that 
they almost cancel. As | x | -* oo the cancellation becomes exact. For 
these reasons it is convenient to divide the range of integration in 
equation (27) into at least four regions : 



< | x | < Z, , 

x t < | x I < V, 

T < \ x\ < X 2 , 

Xa < I x I < oo , 



small | a; |, 
intermediate | x |, 
intermediate | x |, 
large | re |. 



(36) 



From equations (25) , (30) , and (32) we can show that* 



♦These expressions for Q R and Qi may be obtained from the Taylor series 
expansion of the function [sin (x + jy,)]/(x + jy,). 



and 



where 
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Q*(x, y.) = N 2 £ (-!)<-%, J^ , (37) 



Q,(x, y.) =N 2 ± (-l) l+1 A 2(+l ^^ , (38) 



A = s -^ , (39) 

V. 

= cosh y. _ sinhj^ = X (40) 

2/. 2/. iv 

emhfc.^ist, f-1,2.3,..., (41) 

2/. 2/. 



and 



1 


y: l+i 


2£ + 2/c 


+ 1 (2£ + 1)! 


00 


i yV 



A 9M = ^J/--(2*+l) — , 1= 1,2,3, •- . (42) 
It can also be shown that 

*-.- t oTTE 5S z , T^rrn ■ *-i.a.8,-. (43) 
and 

*■ - g ff + a + 1 (§0! ■ * = 1 - 2 ' 3 .-- ^ 

Since the spectrum is an even function of / we can assume without 
loss of generality that 

y. ^ 0. (45) 

For y. ^ 0, from equations (43) and (44) all A t 's are monotonically 
increasing functions of y, , and we can further show that 

< A m+ „ < A ak , k = 0,1,2, ••• , (46) 

and 

^ A„ +I ^ A 2 ,_, , £ = 1,2,3,---- . (47) 

For large y a (for large f/W) , it can also be proved that 

A„ Pd A,-, M S|M. (48) 
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Since it appears that the main contribution to the integral I comes 
from the region of small | x |, assume that | Xi | is small and that 

I = /i + Ib , (49) 

where 

/, = r G(x, y.) dx, (50) 

•'— Xl 

and 

I B = f " G(x, y.) dx + r G(x, y.) dx. (51) 

■» — 00 »Xl 

For small \x\, we have from equations (37) and (38) 

Q R (x)ttN 2 hA 2 x 2 , (52) 

and 

QtixjttN^AsX 3 . (53) 

Let us choose* x t so that exp {-\N 2 A 2 x 2 ) falls to exp (-5) re 0.0067 
for a; = a - ! ,f or that 

10 



.r, = 



ATM. 



(54) 



Since it can be shown from equations (39) -(44) that the minimum 
value of Ao is ^ (at y s = 0) , 

, * (ff. 

Assume that 

x 1 ^ w, (56) 

so that exp { — Q R (x, y,) } is a monotonically decreasing function of x for 
^ x ^ Xi . Equation (56) will be satisfied for all y, if 

N 2 ^ ^ re 3.039. (57) 

Since we are primarily interested in the high-index case, this is not 
a significant restriction. 



* See equation (34) . 

t While there is some degree of arbitrariness in such choices, the bounds 
obtained are not too sensitive to small variations in the value of Xi, x a , and so on. 
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Noticing that 



sin x ^ x, ^ x < 



(58) 



(sin 2 x + sin h 2 y,)* . ( x 2 + sinh 2 ?/. V ,_ Q -. 

(? + !/;)' s I ITT?. I ' (59) 



sin (x + jy.) 
s + jy. 

Since (sinh y,)/(y.) ^ 1, it can be proved from equation (59) that 



sin (x + jy.) 
x + jy. 



< 



sinh ?/. 

y. 



We can show from equations (29) , (50) , and (60) that 
| /, | < 2 f {exp [-Q a (x, y.)] + H(x, y.)} dx, 

Jo 



(60) 



(61) 



where 



[ _ N > *jfcjuji + N > s i^] . (62) 



(63) 



H(x, y.) = exp 
From equation (37) 

From equations (46) and (63) it can be shown that for all y 8 

Q*(x, y.) ^ iV 2 [| A 2 x 2 - ± A 4 * 4 ] , =g x =g x t < (56)*. (64) 
We then have 

| h | ^ 2 f exp {-A^ 2 [| ^ - ^ A 4 z 4 ]} + 2 JH ff(x, y.) Ac. 

(65) 
One can show that 

e 1 ^ 1 + ^-=^ /. ^ / ^ B. (66) 

Since we have 



n « Ar " a « <r 25 A < * 
° = 24^ ="6"XiV% 



^ x ^ x, , 



(67) 
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fOJ? A . 1 1 



f ["25 A 4 1 1 



a; 4 rf.r 



+ 2 f ' tf(a;, 7/.) da; <2 [ exp [-JV 2 £A 2 z 2 ] 

Jo *'o 



exp 



+ 2 



|~25 A 4 1 1 
l_6 A 2 N 2 A 2 ] 



- 1 



Xi 



f° x 4 exp [-N 2 hA 2 x 2 ] dx 
Jo 



+ 2 [ H(x, y.) dx. 

Jo 

From equation (68) it can be shown that 

2/ 



A I < 



#M £ 



(68) 



(69) 



where 



_ 3 / [25 A 4 1 1 

~/5Y I t,t2 sinh w, 

+ 2(-J exp 1 -N -j-* 



1 +iV 



2 sinh y, 
2/. J 



(70) 



Since we know that 

— | p | ^ Re p ^ | p |, V any arbitrary complex number, (71) 

equation (69) gives an upper bound to Re I lm Let us now find a 
lower bound to Re h . 

From equations (29) and (50) we have 



Re 7, = 2 / exp [-Q R (x, y.)} cos [Q^x, y.)] dx 

Jo 

- 2 Re /" exp [-If *£* + jktjl + N 



2 sin (x + jy,) ~\ 
x + jy, J 



dx. 

(72) 



As shown earlier in this paper 



2 Re 



£' GXP [ 



-tf 



gAjh + ^Ji + yin (» + »,) ' | ,,,. 



7 JL * + »• - 



I , r2 sinhw, II, . ,, 2 sinh w, 

^ 2a;, exp -# 2 - ' [\ 1 + A " 



?/. 



?/» 



(73) 
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One can also show that for z real 

z 2 



cos z ^ 1 — 



2 ' 



— oo < z < oo . 



Using equation (74) we can write 

2 / exp [—Q R (x ,%y,)] cos [Qi(x, y,)] dx 

Jo 

> 2 JH exp [-&(*, j,.)][l - ®^)J rfx . 
Now from equations (37) , and (38) we have 

Q B (x, y.) = N 2 [l A 2 x 2 - ± A 4 x{l - U 4 s ^ 



4! 



6! A 4 



±-Ax*(l-^-^A 1 



and 



Qi(x, y.) = iV 2 [| Ay - || A 5 x s -; I - »■ - p •<' 



51 A, 

7! A, 



Q!*"*! 1 1114. */ J 

From equations (76), and (77) one can show that 

Q R (x, y.) ^ hN 2 A 2 x 2 , :g x ^ x, < (30)*, 
and 

Qi(x, y.) ^ \N 2 A 2 x\ ^x ^x, < (42)*. 
Equations (75) , (78) , and (79) yield 

2 / exp [— Q R (x, y.)] cos [Q t (x, y.)] dx 

Jo 

> 2 p exp [-\N 2 A 2 x 2 ] 1 - ^ A 2 x 6 J dx 
= 2 f exp [-AJV 2 A 2 x 2 ] dx-2 p exp [-iJV 2 A 2 x 2 ] dx 



36 A ' 



/* x 6 exp [-£JV 2 A 2 x 2 ]dx. 

•'0 
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(74) 



(75) 



(76) 



(77) 

(78) 

(79) 



(80) 
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Notice that 

f" x* exp [-\N 2 A 2 x 2 ] dx < f x« exp [-%N 2 A 2 x 2 ] dx 

i / ar V 15 (81) 



2\N 2 AJ (N'A 2 y 
We can also show that 

J" exp (- P Y) at < ^i^ 1 . »V H 0. (82) 

We can, therefore, write 

2 J" exp \-\N 2 A 2 x 2 \ dx < (-^J-J e ' 5 - ( 83 ) 

From equations (72)-(73), (75), (80)-(83) it can now be shown 

that 



Re/, >(^):i -^i. (S1 ' 



where 



•exp [-^ Sj ^][l + iV*^} (85) 

We shall now obtain upper and lower bounds to Re I R in equation 
(51) . According to equation (36) let 

I*=h + It (86) 

where 

I 2 = [ " G(x, y.) dx + [ G(x, y.) dx, (87) 

and 

I T = f ' G(x, y.) dx + r ff(s, y.) dr. (88) 

From equation (87) we have 

h\^2 f \ G(x, y.) | dx. (89) 
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Now it can be shown from equations (29) and (60) that 
| G{x, y.) | ^ exp [-Q R {x, y,)] 

From equations (89) and (90) we can write 

\h | ^ 2 f exp [-Q s (x, y.)\ dx 

i «/ \ I »r2 sinli v/,l|, . , r ., sinh y, 

+ 2(tt - x.) exp -A' •" 1 + #" ^ 

L 3/. JL y. 
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From equation (64) 



Qr(x 
where 



, v.) ^ n[1 



A 2 x- - - A <x* 



= N 2 v, Xi ^ x ^ x 



.-^-&^. 



(90) 



(91) 



(92) 



(93) 



It can be shown (see Fig. 8) that v and dv/dx are positive for Xi £ x ^ 
x m = (6) i (A 2 /A i ) i ^ (6)* and that dv/dx is a monotonically increasing 
function of x for .t, ^ £ ^ #„ where 



*- *(£)** *■ 



(94) 



We also know that Q fl (a;, ?/.) is a monotonically increasing function of a; 
for Xi ^ x ^ ir. Let us now assume that x x < V2- 




tfn 



Fig. 8 — Functions v and -g appearing in equation (95). 
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Equations (92) -(94) therefore yield 

2 / exp i—Q R (x, y,)] dx 

•'ii 

< 2 J" exp [-N 2 v] ~dv + 2(tt - V2") exp [-Q R (V2, y.)], (95) 



where 



y, = 


5 . 5 A 4 
' N 2 L 6 A 2 tf : 







2 A 2 \ ' 



(96) 
(97) 



and 

^v, Sv 2 . (98) 

Since we know that 

< §x = j_ i = — iCa (99) 

dv dv ' ' ' 

we can write 



[lL"+-lt^]- 



2 f ' exp [-Aft] * d. < T g?S ; n P exp (-AT 2 *) d. 

J,, ^ l dv J 5i 4 1 •/., 

5 L X 3 A 2 iV 2 A 2 J 

< -= ^ = /"" exp (-iV 2 v) dw. (100) 

r 1 _ 5A 4 1 •>., 
°L 3 A 2 iVSAj 

Since it can be shown that 



/" exp l-p"t] dt = eXp [ -/ V] , p V 0, (101) 



ex p [— p a ] 

2 

from equations (91), (95) and (100), we have 



z\ | < ( Jx)^ 2 ' (102) 
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where 



E 9 = 



(5x)* _ 5 A* _1_ 
3 A, JV 2 A 



exp 



-5(1-14^ l 



6 A, N 2 A. 



+ 2(tt - VS)(^Y exp [-Q,(V5, ?/,)] 



2tt 



+ 2(x - 0,-of^Y exp 



-JV 



2 sinh 2/« 



Z/. - 



1 + JV 



, sinh y 



Xx < V2". 



■I 



(103) 



Similarly, if .?! > s/2, one can show that 
E 2 = 2(r - zof^T^Vlexp [-Q,(x, , 2/,)] 



2t / 
+ exp 



^j 



_^ 2 sjnh ll || 1+iv2 sinh^ 



(104) 



y. JL 2/. _ 

Let us now consider the range of integration t < x < » . For z ]» ?/. 

exp [-#*(&, 2/.) + jQi(x, y.)] 



[" /-..osinh ?/» , r2 i ^ '" •»' 
exp - ( JV" - — — N cosh y 



.'/a 



+ j[\x + AT sinh y. 



COS X 

X 



(105) 



and 
exp 



[-N-'-^ + iK*] 



2 sin (a; + jy 
x + ;V 



[ 



1 + JV 



M — [ 



. t2 sinh ?y, . A 



2/. 



■] 



, . -.alsin x , . . cos a; . , 
1 + JV cosh ?/. + 3 smn 2/- 

|_ X x 



(106) 



Let us choose the point x 2 + jV« along the path of integration so 
that the amplitudes of the two terms in equations (105) and (106) 
differ by less than 10.5 percent [(JV 2 cosh y,)/x 2 ^ 0.1]and their relative 
angle departs from 180° by less than 0.1 radian [(JV 2 sinh y.)/x 2 ^ 0.1]. 
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Such a point x 2 + jy. is given by 

x 2 = ION 2 cosh y, . 

We assume that x 2 > t ^ Zj , or that 



N 2 > fQ^ 0.31416. 



(107) 



(108) 



Since from equation (57) N 2 ^ 30/V 2 , tliis inequality is always satisfied. 
We shall now write 



It = h + h, 



where 



and 



h = J G(x, y.) dx + f G(x, y.) dx, 
h= I ' G(x, y.) dx 4- [ G(x, y.) dx. 



Noticing that 

I sin (x + jy.) \ < cosh y, < cosh y, 
I x + jy. | = (x 2 + y 2 )* - x 

we can show that 



(109) 



(HO) 



(HI) 



(112) 



7 3 | ^ 2 f exp [-Q R (x, y.)] dx 

t.t-2 sinh y, \ 
+ 2(x 2 -r) exp [-AT — ^J 



In 


(-)' 




W / 


x 2 


— 7T_ 



1 + N 2 cosh 7/. — — , (113) 



and 



exp [-Q*(3, y.)] dx ^ £ / exp [-Q R (z, ?/,)] cfo 

-r E f 'exp [-Q«(rc, y.)] dx, (114) 

*=1 •'2/r 



where K is an integer such that 

(2K + 1)tt > x 2 ^ (2K - 1)tt. 
One can show that equation (115) is satisfied if 



(115) 



SPECTRAL DENSITY BOUNDS 



793 



A' = INT 



BL + I 

2t T 2 



or 



g -/yr[ ^^ + |]- (ii6) 



We now have 
(wir) a 



< 



s < 



< 



(n + 1)V 



M 5 + tf = * 8 + v\ = (n + 1)V + y\ ' 

jitt ^ X £ (n + !)«-, n = 1, 2,3, 



(117) 



and 

1 



y, sin a; 2 1 — cos .r . 
tanh ?/. .t x 2 ~~ 



{21 - 1)tt ^ X ^ 2/tt, ^ = 1, 2, 3, • • ■ . 
From equations (31), (117), and (118) we can now prove that 



(118) 



2 / exp [— Q a (x, ?/.)] dx 

{-l J 121-1) x 



K 

I 

l-\ 



^ 23 ir exp 



-N 



2 sinh ?/„ 



1 



y * 1+-*-* 



[ xr2 sinh y.l yr^ 
— N* — 2s exp 
2/. J <-i 



7T (2* - 1) _ 

N 2 y, sinh ?/, 



(2( - if + h 



(119) 



Further it can be shown that* 



A [N 2 y. sinh «. 1 1 

g exp |_ ? ___-^J 

< exp r ^-aiDhy, 1_| 



+ (K - 1) exp 



iV 2 ?/, sinh y, 



9 + 



rf/»\r 



(120) 



*The upper bound derived in equation (120) can be improved in various 
ways. Since this makes only a minor contribution to the total integral we shall 
be satisfied with this simple bound. 
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From equations (119) and (120) we have 
XI / > exp [-Q K (x, y.)] dx 



< 



it exp 



-N'- 



sinh y, 
y. J 



exp 



[ N 2 y, sinh y, 1 j 

1 + (2/.A) 2 J 



■ ITT in \ N 2 y t sinh y, 1 

+ (tf-l)exp[_ 7 ^^ 



) 2 . 
For 2ftr + tt/2 g x ^ {21 + 1)tt, / ^ 1, we can show that 

?/, sin re , 2 1 — cos x 



tanh y, x 



> 1 - 



?/, 



tanh ?/, (4£ + l)ir 
In this range of a; 



sin x + 



y. 



(2/ + 1)V 



2 (1 — cos re) 



cos re rjs 0, 
sin re ^ 0, 



and 



dJ e 



y. 



cos x + 



?y; 



da; tanh y, (U + 1)tt 

We, therefore, have 

1 — cos re 



{2t + 1)V 



y, sin re , 
tanh ?/, re 



where 



VM = 1 - 



y. 



+ 



tanh ?/. (4/ + l)ir (2£ + 1) V 



It can be shown that 



(121) 



(122) 

(123) 
(124) 

sin re ^ 0. (125) 



^ J t (x) ^ J t \2hc + |) m V t (y t ), (126) 



(127) 



V<(0) = 1 - 



(U + 1J7T ' 



(128) 
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that V ( {y,) reaches its minimum at 

»■ " [ l + «ti]" (129) 

and 

[F«(V.)U « 1 - | [l + ^r7_j" > 0, 12 1. (130) 
Now for 2/V ^ .r. ^ 2/V + x/2, it can also be shown that 

i y> sm 31 i 2 1 ~ cos -t y, sin .r 

tanh ?/, a; re 3 tanh y, 2lir 

+ (4/ ffi )V a (1 - cos .r) - £,(*). (131) 
One can prove that L ( (x) reaches its minimum at 



x = 26r + tan" 
and 



(41 + 1) 
. 8C y 



■ \ I , (132) 

, tanh ?/, J 



[LXaOU - CT«fo.) - i + (4g + l)v 



l_(4* + 1 1)V 4 + «V tanh^7/J (133) 
Next we can show that 

C/,(0) - 1 - JL < 1 - ^^ = 7,(0), (134) 

that U ( (y.) is a monotonically decreasing function of y, (see Fig. 9), and 

Lim U t (y.) = 1 - (4 ' 3 y > 0, * > 1. (135) 

It can aleo be proved by numerical methods (see Fig. 9) that 

U<ty.) ^ UM, y.^0, f^l, (136) 
and 

V t (y.) > UM.y. £ 0, I ^ 1. (137) 
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We therefore conclude that 

« y, sin x . 2 1 — cos x TT . N . . 4y;' 

L62&r 4 + 4r'tanh 2 tfJ' " *' ' ,8U ' 

2tw ^ x ^ (2C+ l)ir. (138) 
Equations (31), (117), and (138) show that 

. . 2 sinh y. rT , . (2QV 

<?,(*, ?/.) ^ # -^— I7.df0 (20V + y * . 

26r g .r £ (2£ + 1)tt. (139) 
From equation (139) we can now write 

2 / exp [-Qb(x, ?/,)] da: 

^ v- T ^ aiahy. „, . (2QV 1 

= Texp [_^5^ a(!/i) j 

It can be shown that 

A [ N'y, sinh y. r , , , 1 1 

+ (X _ 1) ^[^^^. ) _^]. (141) 
Equations (140) and (141) yield 

2 / exp [-Q R (.r, y.)] <fc < tt exp -iV 2 S -^ £/,(</.) 

<-l •'2<T L </« J 

+ (K _ 1) exp [> K^jUb F , w ^L^]}. (142) 
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From equations (113), (114), (121), and (142) we can write 



/» I < 



2tt 



N 2 A. 



E 3 , 



where 



E> = 



.y^ „[_,-!»] 



r (N 2 y, sinh v, 1 



i V 2 A 

2tt 



+ (X - 1) exp 



iV 2 ?/, sinh y. 



1 



9 + (y.A)V_ 



j"\ r2 y, sinhy 



fc.A)*J 



4- (tf - 1) exp N 



4 + (2/ 

ffi(v.) 



Iu|S 



.To — 7T 



If 2 is a complex variable, it can be shown that 



exp (2) — 1 — 2 [ ^ 



2 



exp \z \. 



From equations (112), (145) -(147), we can write 



(143) 



16 + (y.A) 1 . . 

+ 2(x 2 - tt) exp r-AT 2S -^jM 1 + N 2 cosh?/, 
Finally, from equation (111) we have 

I U I ^ 2 j" I G(x, y.) I <fe. 
Now from equations (27) -(28) 

|^ l2/ .)| = exp[-^ 2S -^] 

J r^ 8m(x + ^) 1 _ [\ + N2 sm(x + jy.) -\ I 

I L a; + j?/. J L a; + jy. J I 



(144) 



(145) 



(147) 
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I /« I £ exp r-JV 2S -^^V cosh 2 y. 



dx 



r 1 l~N 2 cosh y.~] 

N* cosh 2 y, f , r2 sinh y.~] 
= — exp —N 

y. L y. J 

•>0 



N* cosh 2 i/. 



exp 






2/. 



N 4 cosh 2 y, _ 2 V 

y 2 > 



r^sinhy.l 

V N 2 cosh y, _. 

L y. 



(// 



exp 



sin 



d0. (148) 



Since 



^ sinifl ^ 0, 0^0^ sin" 
we can show from equation (148) that 



Mi ^ Z 
{xl + 2/ 2 )' < 2 



T - „ A 7 " 4 cosh 2 y, 
I* I < *- exp 



2/. 



-AT 



sinh y, 



2/. 



./;""' •■'"■'*•"" exp [^^e]de 



= A 7 " 2 cosh y, exp 



["_ ^Binhy.1 



• jexp I 



y. 



, r2 cosh y, . _i 

A 7 •" sin , 2 2 

2/. (.x-2 + I/.) 



,] - .}■ 



Now we have 



^ sin" 1 cr g|<r, 0i£ <r £ 1. 



From equations (149) -(151) , we can write 
| 7 4 | < N 2 cosh y. exp [-iV 2 s -^ h -^] 



(149) 



(150) 



(151) 
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where 
E, = 



A' 2 cosh y, exp [-i\T S -^^j[exp (£) - l]- (153) 



2tt 



From equations (27), (49), (69), (71), (84), (86), (102), (109), 
(143), and (152) we can write the following bounds for Re I: 



(j£ju -«-*-*-*! 



2tt 



< Re / < (j^-J [1+E l +E a +E z + E 4 ]. 

It has been shown that 

sinh y, _ 2 / 
y. y. N*W' 

V. UPPER AND LOWER BOUNDS TO S v {f) 

We have shown in the previous section that 

Sy(fi = exp (-W*(/) + £= [u-.tf + IF) - u. x (j - W)} 



2 IP 



+ 



^expl^L^^- 8 -^ 



2rW 



y, - 



where 



Gra-jv -«-*-*-*! 



< m < (j^j-J [1 + -E\ + # 2 + E a + *J, 

cosh ?/, sinh y, _ f 
~y. y, ~ N 2 W ' 

_ sinh y, _ 2 / 
2 " 2/. W.iV'TP 



Parameter 
Ei 
E{ 
E, 
Ex 
E 4 



Equation 
70 

85 

103 or 104 
144 
153 



(154) 
(155) 



H (156) 



(157) 

(25) 

(155) 
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For # 2 = 10 and 25 we plot, in Figs. 10-16, E l , E[ ,E 2 ,E 3 , and E 4 . 
Notice that E 1 , E[ , E 2 , E a , and E t appearing in these bounds are 
all very small compared to unity so long as the modulation index is 
moderately high, and that E x ,E Z , and E A are monotonically decreasing 
functions of / and N 2 . Also notice that E[ and E 2 may first increase 
(see Figs. 11, 12) with y. (or /), reach their maxima and then decrease 
with y, * It can be shown that these maxima are all very small com- 
pared to unity for all N 2 which are even moderately high. 

For all /, we can then write 

fe)^ 1 - *> < ' < fe)' (1 + D) ' (158) 



where 



and 



C = E[ + E 2 + E 3 + E A , (159) 



D = E l + E 2 + E* + E A . (160) 

From Figs. 10-16 and expressions for C and D, we can show that 
C and D are both small (<2%) compared to unity for N 2 > 25 and 
for all /. Hence we deduce that 

: - 7 ' l (161) 



kN 2 A_ 

and that the fractional error in this approximation is very much 
less than unity ( <2%) . 

For N 2 — 10 and 25 the spectral density S v (f) and the fractional 
error C and D obtained from equations (158) -(161) are plotted in 
Figs. 17-20. From these figures notice that C and D are less than 
10 percent for iV 2 > 10,f and that 

C < 2%, for JV 2 ^ 25, (162) 

D < 2%, for JV 2 ^ 25, (163) 

proving the assertion made earlier in this paper. 

For N 2 = 6 the spectral density obtained from equations (158), and 
(161) is given in Fig. 21; the percentage error between this spectral 



* One of the terms in Ex is independent of / and N 2 . 

t By modifying the contour of integration we also have been able to show 
that C and D are less than 8% for N" ^ 10. Since this modified contour leads to 
unnecessary complications, we have not given that modified contour analysis in 
this paper. 
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Fig. 9 — Functions U t (y,) and F,(i/.1. It can be observed that F,(j/ 3 ) > C/i(t/») > 0, 
I 2 1. 
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Fig. 10 — Parameter E t as a function of hjyi w" 
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Fig. 11 — Parameter Ei as a function of I 2jp\ W 
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Fig. 12 — Parameter £2 as a function of o/V* W 
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Fig. 13 — Parameter E% as a function of oA^ W> w ^ n ^ 2 = 10. 
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density and that obtained from equation (19) has been plotted in 
Fig. 22 (for a set of values of f/W) . The scatter diagram in Fig. 22 
indicates that the spectral densities obtained from the two methods 
agree veiy closely and that the saddle-point approximation error is 
not related in a simple way to the truncation error (it does not seem 
possible to draw a smooth curve through the points shown in Fig. 22) . 

For all practical purposes, including interference calculations, esti- 
mation of the spectrum to such an accuracy is almost always suf- 
ficient. It can therefore be said that the saddle-point approximation 
given by equations (25), (155), (158), and (161) is a good approxima- 
tion to S v (f) as long as the modulation index is even moderately high 
(N 2 > 10). The spectrum can be estimated by this method for all 
values of / even when it is millions of decibels smaller than the con- 
tinuous part of the spectrum at / = 0. 

Now compare the spectrum obtained from the quasistatic approxi- 
mation* to that obtained from saddle-point approximation. For this 
purpose, the spectra obtained from equation (21) for N 2 = 10 and 25 are 
plotted in Figs. 17, and 19. We see that the spectra obtained from the 
quasistatic approximation agree very closely with those obtained from 
the saddle-point approximation for low frequencies, but that the quasi- 



* See equation (21). 
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Fig. 17 — Spectral density of an angle-modulated wave, with gaussian phase 
modulation with a rectangular spectrum. N = (10) 1 " zz 3.162 radians, rms phase 
deviation. 



static approximation to S v (j) is too small for large /.* In fact for N z = 10 
the quasistatic approximation is 30 dB too small for f/W « 13.5. We 
have therefore shown that the quasistatic approximation to the spec- 
trum cannot be used in any interference calculations or in any calcula- 
tions where the behavior of the spectrum on the tails is of importance.! 
The saddle-point approximation can be used at all frequencies as long 
as N 2 is moderately high. 

* For small / (or small y.) it can easily be shown that the saddle-point ap- 
proximation reduces to the quasistatic approximation. 

t The higher the rms phase deviation, the further out will the low-frequency 
(quasistatic) approximation be valid, 
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Fig. 18 — Bounds on fractional error in saddle-point approximation to the 
spectrum. N 2 = 10. 



VI. RESULTS AND CONCLUSIONS 

A simple method (called the saddle-point method) has been pre- 
sented in this paper to estimate the spectrum of a sinusoidal carrier 
phase modulated by gaussian noise having a rectangular power 
spectrum. 

This method gives upper and lower bounds to the spectrum and 
shows that these bounds are very close for all / and for all moderately 
high phase deviations. We also show that the fractional error in the 
saddle-point approximation is less than 2 percent for N 2 ^ 25 and 
for all /. 

The calculation of the spectrum by the saddle-point method is 
rather simple. For a given value of /, N 2 , and W, we calculate y a 
from equation (25) and A 2 from equation (155). The spectrum S v (f) 
is then calculated from equations (156) and (161). 

We have also shown in this paper that the quasistatic approxima- 
tion to Svtf) is only good at low frequencies, and that for large / the 
results obtained from that approximation are too small. 
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APPENDIX 



It can be shown (see Ref. 7, p. 114) that 

if" T n 2 r w 1 

S v (f) = j^pp J ex P ~2W j_ w (- 1 - cos 2t ^ t ) d* 1 

•exp [— j2ttJt] dr, 



(164) 
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Fig. 19 — Spectral density of an angle-modulated wave, with gaussian phase 
modulation with a rectangular spectrum. N = 5 radians, rms phase deviation. 
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Fig. 20 — Bounds on fractional error in saddle-point approximation to the 
spectrum. N 1 = 25. 
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Fig. 21 — Spectral density of an angle-modulated wave, with gaussian phase 
modulation with a rectangular spectrum. N = (6) 1/2 SS 2.449 radians, rms phase 
deviation. 
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Fig. 22 — Percentage error between the spectral densities obtained from saddle- 
point approximation and that obtained from equation (19). It does not seem 
possible to draw a smooth curve through the points shown in this figure. We 
have, therefore, shown the error as a scatter diagram. 



or 



^-/-".{^["^/-K 1 "" "»#')*] 



- exp [-TV 2 ] 1 + :\ 

From equation (165) , it can be shown that 
. I . 2 sinhj/."| 

, = exp L-* _j.j 



, sinp I jxp 



V 



■e' sv dp. (165) 



, sin (x + jy.) 



cos ^ (a: + J2/ 



0*»] 



- 1 + N* 



x + jV, . 



e' Xl rfx, 



(166) 



and 



Q.(b, V.) = N 2 S -^p - Re |£ /_" cos ^ (z + i?/.) d Ml (167) 
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or 

Q R (x, y.) = N 2 S -^ - |p [ cos £ x cosh £ y. d M . (168) 



y 

Equation (168) yields 



dQ*(x. y.) _ tf { w JL^JL 
dx " W 



I w sm w xco ^w y ' dfl - (169) 



For ^ vx/W ^ *-, 



sin^x^O. (170) 



For y, ^ 0, and ^ n ^ IF, 

^ a; ^ 7T, (171) 



and from (169), 



W s ' m w xcostl w y ' = > (172) 



^# y ' } ^ 0, ^ x g tt. (173) 



From equation (173) we then conclude that Q R (x, y.) is a mono- 
tonically increasing function of x for ^ x ^ ir. 
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